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C ~ ~) • Abstract. 

Faint extragalactic radio sources provide important information about the global 
C ~ ~) ' history of star formation. Sensitive radio observations of the Hubble Deep Field 

, and other fields have found that sub-mjy radio sources are predominantly associated 

■ with star formation activity rather than AGN. Radio observations of star forming 

galaxies have the advantage of being independent of extinction by dust. We use the 
FIR-radio correlation to compare the radio and FIR backgrounds, and make several 
conclusions about the star forming galaxies producing the FIR background. We then 
use the redshift distribution of faint radio sources to determine the evolution of the 
radio luminosity function, and thus estimate the star formation density as a function 
of redshift. 



1. Introduction 

Faint radio sources provide important information about global star formation history. 
Sensitive radio observations of the Hubble Deep Field (HDF) (Richards et al 1998, 
0) and other fields well-studied at optical wavelengths (Windhorst et al 1995, Q, 
Fomalont et al 1991, Q) have shown that sub-mjy radio sources are predominantly 
associated with star formation activity rather than active galactic nuclei (AGN). The 
radio luminosity of a galaxy is a reliable predictor of the star formation rate (SFR) 
for local galaxies (Condon 1992, Q Cram et al 1998, ||). Estimates of star formation 
based on radio observations also have the advantage of being independent of extinction 
by dust, which has caused much difficulty in the determination of star formation 
history from optical data. 

In section 2, we make use of the tight correlation between radio and FIR luminosity 
for star forming galaxies to compare the FIR and radio backgrounds and to study the 
sources producing both. In section 3, we determine the evolving radio luminosity 
function from the observed redshift distribution of faint radio sources, and then 
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estimate the history of star formation to a redshift of about 3. Throughout, we 
assume f2 TO = 1, f^A = 0, and 7Jo=50kms _1 Mpc _1 . We define the radio spectral 
index as S u oc v~ a . 

2. FIR vs. Radio Backgrounds 

This section follows our recent paper (Haarsma & Partridge 1998, |J). The far infrared 
(FIR) background was recently detected with DIRBE (Hauser et al 1998, fjj; Dwek et 
al 1998, [|J), and is most likely the collective emission of star forming galaxies. We use 
the radio-FIR correlation for individual galaxies (Helou, Soifer, & Rowan-Robinson 
1985, H) to calculate the radio background associated with the FIR background, 
assuming that the bulk of emission is from z ~ 1. We find the radio background 
associated with the FIR background has a brightness temperature of T 4 ocm = 0.31 K, 
or T^ocm ~ 15 K (scaled using a spectral index of a — 0.7). At 170cm (178 MHz), 
the observed radio background is T 170cm = 30 ± 7 K (Bridle 1967, |h}). This allows 
us to draw several conclusions about the faint sources making up the FIR background: 

(i) The radio emission from these sources makes up about half of the observed 
extragalactic radio background. (The other half is the summed radio emission of 
AGN.) 

(ii) Since (i) is in agreement with other radio observations (Condon 1989, [0), the 
FIR-radio correlation appears to hold even for the very faint sources making 
up the FIR background. This confirms the assumption that the FIR background 
between about 140 and 240 /im is dominated by star-formation, not AGN activity. 

(iii) By quantitatively comparing the radio and FIR backgrounds, we find a 
relationship for the sources contributing to the background, 

a 0^y) =°- 20±0 - 05 ' w 

where a is their radio spectral index, A is the fraction of the radio background 
they produce (from (i), A ~ 0.5), and z is their mean redshift. This function is 
plotted in Figure [|. Note that the redshift z is the mean redshift of the sources 
dominating the FIR and radio backgrounds, which is not necessarily the redshift 
of peak star formation activity (see §3). 

(iv) By extrapolating the 3.6 cm logiV — log S curve to fainter flux densities, we 
estimate that most of the FIR background is produced by sources whose 3.6 cm 
flux density is greater than about 1 /iJy. This lower limit is consistent with 
other work (Windhorst et al 1993, |l2|), but has more interesting observational 
consequences. An RMS sensitivity of 1.5 /iJy has already been reached in VLA 
observations (Partridge et al 1997, |r|). The log N — logS curve indicates that 
the number density of S > 1 /iJy sources is about 25/arcmin 2 , similar to some 
model predictions (Guiderdoni et al 1998, At this density, these sources 
will cause SIRTF to encounter confusion problems at 160 fj.ui. 



3 



3. Radio Star Formation History 

In this section, we use the redshift distribution of faint radio sources to determine the 
evolution of the radio luminosity function, and the evolution of the star formation rate 
density. 

3.1. Data 

Three fields have been observed to microJy sensitivity at centimeter wavelengths 
and also have extensive photometric and spectroscopic data: the Hubble Deep Field 
(HDF), the Medium Deep Survey (MDS), and the V15 field. Table [l] gives the details 
of the three fields and references. For the first time we have a sample of micro Jy 
radio sources with nearly complete optical identifications and about 50% complete 
redshift measurements. We assume that all sources detected at these flux levels are 
star-forming galaxies, since optical identifications indicate that ^80% of these radio 
sources have spiral or irregular counterparts (Richards et al 1998, Q). The known 
quasars (two in the MDS sample, none in the HDF or V15 samples) were removed. In 
the flanking fields of the HDF, we have used the relationship between redshift and K- 
band magnitude (Lilly, Longair, & Allington-Smith 1985, ||l5| ) to estimate redshifts for 
10 sources without spectroscopic values. For the remaining sources without redshifts, 
we arbitrarily selected redshifts to fill in gaps in the redshift distribution, in order 
to illustrate the total number of sources that will ultimately appear on the plot. 
Photometric redshifts for these sources are currently being calculated (Waddington & 
Windhorst, in preparation), and will be included in future work. 

To compare these data to the model, we calculate n(z), the average number of 
sources per arcmin 2 in each redshift bin. This requires a correction for the varying 
sensitivity across the primary beam of the radio observations (Katgert, Oort, & 
Windhorst 1988, @; Martin, Partridge, & Rood 1980, 0). For example, a faint 
source which could only be detected at the center of the field contributes more to n(z) 
than a strong source which could be detected over the entire primary beam area. The 
resulting redshift distributions are plotted in Figures [| and ||. 

It is interesting how different the MDS and HDF distributions are, even though 
the surveys were both performed at 8 GHz with similar flux limits. The average 
source density (including all sources) is 1.26 sources/ arcmin 2 in the HDF, but 2.63 
sources/arcmin 2 in the MDS field (the V15 field at 5 GHz has 0.736 sources/arcmin 2 ). 
The density of sources in the MDS field is over twice that of the HDF field, possibly 
due to galaxy clustering. In the analysis below, we fit the model to the data in all 
three fields simultaneously. 

3.2. Calculations 

In order to determine the star formation history, we must first determine the evolving 
radio luminosity function. We used two versions of the local 1.4 GHz luminosity 
function for star-forming/spiral galaxies. We define the luminosity function <fi(L e ,iA) 
as the number per comoving Mpc 3 per d log 10 L °f star- forming radio sources with 
emitted luminosity L e ,i.4(W/Hz)at 1.4 GHz. Condon (1989, pjj|) uses the following 
form for the luminosity function (but different notation), 

log 10 [<£(Le,i.4)]dlog 10 L = 28.43 + Y - 1.5 log 10 L e ,i.4 
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Field 


Location 


Band 


Flux limit 


N 


N z Reference 


Hubble Deep Field 


12hH 


-62d 


8 GHz 


9/jJy 


29 


13 (+10) H 


Medium Deep Survey 


13hH 


-42d 


8 GHz 


8.8/xJy 


19 




VI 5 field 


14hH 


-52d 


5 GHz 


16 /zJy 


35 


18 10 



Table 1. Summary of three deep radio surveys. The flux limit of the complete 
catalog for each field is approximately 5 times the RMS noise of the observation, but 
varies across the field. N is the total number of sources above the flux limit, and 
N z is the number of those sources with spectroscopic redshifts. An additional 10 
redshifts were estimated in the HDF from their K-band magnitudes. 
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d\og w L, 
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with the fitted parameters for star-forming galaxies of Y = 2.88, X — 22.40, W = 2/3, 
and B — 1.5. Serjeant et al (1998, p8|) use the standard Schechter form, 



cj>{L e ^ A )d\og w L = 4>* In 10 



L 



e,1.4 



(l+m) 



exp 



d\og w L (3) 



where a factor of L In 10 has been included to convert the function from dL to dlog 10 L. 
Serjeant et al find fitted parameters of <j>* = 4.9 x 10~ 4 Mpc -3 , L* = 2.8 x 10 22 W/Hz, 
and ai = —1.29. 

To describe the evolution of the luminosity function, we use the functional form 
suggested by Condon (1984a, |19 , eq. 24), a power-law in (1 + z) with an exponential 
cut-off at high redshift. The luminosity evolves as 



f(z) = (1 + z) Q exp 
and the number density evolves as 
g(z) = (1 + z) p exp 



(4) 



(5) 



This gives six parameters {Q, q, z q , P,p, z p } to use in describing the evolution. When 
fitting for the parameters, we constrained the functions g(z) and f(z) to the physically 
reasonable ranges of 1 < g(z) < 100, 1 < f(z) < 100 for < z < 3. The general 
expression for the evolving luminosity function is then (Condon 1984b, pcj|) 



<t>( L e,iA,z) = g{z) 



L e ,lA 



(G) 



To use this expression at an arbitrary observing frequency v and redshift z, we must 
convert the observed luminosity L a „ to 1.4 GHz and do the K-correction, i.e. 



(7) 



where a is the radio spectral index, as defined in §1. We have assumed a — 0.4 for all 
calculations in §3 (Windhorst et al 1993, JT1). 
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The evolving luminosity function can be used to predict the observed redshift 
distribution. The number of sources per redshift bin Az that could be detected in a 
survey of angular area AQ and flux limit Su m at frequency v is 

/•inf 

n(z) = V c (z, Az, Ail) / <j){L e . 1A ,z)d\og 10 L (8) 
Jl'(z) 

where the lower limit of the integral is 



Dl{zY 2 
Mpc 



(9) 



and Dl is the luminosity distance. The comoving volume in a shell from z to z + Az 
and angular size An is 



V c (z, Az, Afi) = J dri J r 2 {z)dr 

An / ster \ [r 3 (z + Az) - r 3 (z)] 



stcr \ 1.18 x 10 7 arcmin 2 

where the comoving distance is 



(10) 



for our assumed cosmology (see §1). We have used An = 1 arcmin 2 for comparison 
to the data in Figures ^ and ||. 

The evolving luminosity function also allows us to calculate the star formation 
history. For an individual galaxy, the star formation rate is directly proportional to 
its radio luminosity (Condon 1992, Q): 

SFR = q( wj)— (12) 

\5-3 x 10 2 * ( ck)"°- 8 + 5.5 x 10 2 " (c^r V yr 

The radio luminosity is primarily due to synchrotron emission from supernova 
remnants (the first term in the denominator) plus a small thermal component (the 
second term) . Both components are proportional to the formation rate of high-mass 
stars which produce supernova (M > 5M Q ), so the factor Q is included to account 
for the mass of all stars (0.1 — lOOM©), 

n g M\jj{M)dM 

where ip(M) cx M~ x is the initial mass function (IMF). We have assumed throughout 
a Salpeter IMF (x — 2.35), for which Q = 5.5. If an upper limit of 125Af Q is used, 
then Q = 5.9. 



In order to use eq. 12 at high redshift, both L v and v in the equation must be 
K-corrected to the emitted luminosity at the emission frequency. Are there other ways 
in which this relation evolves? The thermal term is much smaller than the synchrotron 
term, so evolution in the thermal term will have little effect. In the synchrotron term, 
the dependence on the supernova environment is weak. One component that might 
cause significant evolution in eq. O is an evolving IMF, entering through the factor Q. 
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In active starbursts, the IMF may be weighted to high-mass stars (Elmegreen 1998, 
pj|), which would result in a smaller value of Q. However, the smallest Q is unity 
(when virtually all mass occurs in high-mass stars), so the strongest decrease due to 
IMF evolution would be roughly a factor of five. 

To determine the star formation rate per comoving volume, we simply substitute 
the radio luminosity density for L v in eq. |l2|. The star formation rate depends on the 
emitted (rather than observed) luminosity density. The luminosity density emitted at 
1.4 GHz can be easily found from the evolving luminosity function, 

/inf 
L e _ 1A (j)(L e:1A ,z)d\og 10 L. (14) 
-inf 

Thus the predicted star formation history is 

where 1.4 GHz is used in the denominator of eq. |lj (no K-correction is needed because 
the luminosity density is the emitted value). 



3.3. Results 

We use the formulation of §3.2 to determine the star formation history from the 
evolving luminosity function. To determine the evolution parameters, we compare 
the model to the observed n(z) for the three surveys. We immediately found that 
pure luminosity evolution [f(z) = (1 + z) 3 and g(z) — 1], as often suggested in the 
literature, is a poor fit for the faint star-forming galaxy population (the predicted 
n(z) is too small and has a very long high-redshift tail). The model fit of Condon 
(1984a, |ll|), {Q = 3.5, P = 1.75, p = 1.8, z p — 1} with no exponential cut off in 
luminosity evolution, is much better (more reasonable redshift dependence, but n(z) 
is still too low). To improve on these models, we adjust the evolution parameters 
{Q, q, z q , P,p, z p } to improve the model fit to the n(z) data, using a downhill simplex 
algorithm (Press et al 1992, |22|]) to find the global x 2 minimum. We performed this 
fit using both the Condon (1989, |]l~l| ) luminosity function (Model C, see eq. ||) and 
the Serjeant et al (1998, Q) luminosity function (Model S, see eq. 

In Model C we use the luminosity function of Condon (1989, |ll|). The fitted 
evolution parameters were {Q = 7.6, q — 1.3, z q = 0.48, P = 1.6, p — 1.2, z p = 1.8}. 
The resulting evolution factors f(z) and g(z) are plotted in Figure || and the resulting 
luminosity function is shown in Figure |3j. Although the term (1 + z) 7 6 seems 
extreme, when combined with the exponential cut-off the luminosity evolution f(z) 
is reasonable. The fit to the the redshift distribution is shown in Figure ||. The 
fit significantly underestimates the total number of sources in the MDS field, but 
only slightly underestimates the other two fields. The V15 survey has the largest total 
number of sources and thus has the most weight during fitting, so the result is a better 
fit for V15 than the other fields. Finally, Figure || shows our predicted star formation 
history (heavy line) along with model predictions from several others (thin lines). The 
vertical lines indicate the 1/VN uncertainty, where N is the sum of galaxies at that 
redshift from the three surveys. The Model C prediction is in good agreement with 
other models at low redshift (the curve follows closely the prediction of Pei & Fall 
1995, Q, as plotted in Dwek et al 1998, ]^], figure 3), which is impressive given that 



7 



no free parameters were adjusted to fit the z = value. The predicted star formation 
history peaks around a redshift of 1, and falls off more quickly than other models at 
high redshift. 

In Model S we use the luminosity function of Serjeant et al (1998, |Q) (see eq. || 
above). The fitted evolution parameters were {Q — 4.3, q = 2.1, z q = 1.3, P = 1.3, p = 
1.7, z p = 2.3}. The resulting evolution factors f(z) and g(z) are plotted in Figure || 
and the resulting luminosity function is shown in Figure |7]. Despite very different 
individual parameters (Q = 4.3 vs. Q = 7.6), the two fits have similar functions f(z) 
and g(z). The predicted redshift distribution (Figure ||) is peaked at a slightly lower 
redshift and has a slightly longer tail then Model C. The predicted star formation 
history (Figure ^|) has a larger local value than Model C, but still less than that 
predicted by Baugh et al (1998, [p4[ ) (thin solid line). The peak is around a redshift 
of 1.4, and falls off less rapidly than Model C at high redshift. 

3.4- Discussion 

The star formation histories predicted by Model C and Model S both fall off 
more quickly at high redshift than model predictions by others. However, we are 
considering several refinements to our model that might modify this result. We are 
currently determining additional photometric redshifts (Waddington & Windhorst, in 
preparation), which will make the modeling more reliable particularly at high redshift. 
The predicted shape of the star formation history is limited by the functional form 
we chose for evolution (eq. ^ and [5]), and we plan to experiment with other functions. 
If the IMF is evolving, or is dependent on environment, this would also affect our 
results. The relationship between star formation rate and radio luminosity (eq. [l2|) 
might be evolving in addition to its dependence on an evolving IMF. Finally, we have 
not explored the dependence of our results on cosmological parameters. 

This method has the potential to be an important indicator of star formation 
history. Radio luminosity is a reliable indicator of star formation rate in local galaxies, 
and is not affected by dust extinction. While others are performing similar calculations 
(Cram et al 1998, [|; Cram 1998, f|; Mobasher et al 1999, |||, Serjeant et al 
1998, [[l8)), the survey data used here are complete to a substantially lower flux limit, 
with nearly complete knowledge of optical counterparts and ~50% completeness in 
redshifts. This allows us to place stronger constraints on the evolving radio luminosity 
function and to probe star formation activity to much higher redshifts. 
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Figure 1. Relationship between the radio spectral index a, the ratio of star- 
formation flux to the total radio background A, and the typical redshift z for the 
sources making up the FIR, background. 
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Figure 2. Evolution functions for Model C. The solid line is f(z) (luminosity 
evolution), and the dashed line is g(z) (number density evolution). 
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Figure 3. Evolving luminosity function for Model C. The labels indicate the redshift 
for each curve, for z = 0, 0.5, 1.5, and 2.5. 
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Figure 4. Redshift distribution n(z), number per arcmin 2 per redshift bin. The 
curve shows the prediction of Model C. The data histograms indicate spectroscopic 
redshifts [heavy line] , spectroscopic redshifts plus estimates from K magnitudes [thin 
line histogram], and all sources (including arbitrary redshifts for remaining sources) 
[dotted line]. 




Figure 5. Star formation history. The heavy curve is the prediction of Model C, with 
vertical lines indicating the uncertainty. The remaining lines indicate star formation 
histories predicted by several other models: solid line (Baugh et al 1998, [22]), dotted 
line (Guiderdoni et al 1999 [26]), other broken lines (Dwek et al 1998, [8], figure 3). 
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Figure 6. Evolution functions for Model S. The solid line is f(z) (luminosity 
evolution), and the dashed line is g(z) (number density evolution). 
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Figure 7. Evolving luminosity function for Model S. The labels indicate the redshift 
for each curve, for 2 = 0, 0.5, 1.5, and 2.5. 
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Figure 8. Redshift distribution n(z), number per arcmin 2 per redshift bin. The 
curve shows the prediction of Model S. The data histograms indicate spectroscopic 
redshifts [heavy line] , spectroscopic redshifts plus estimates from K magnitudes [thin 
line histogram], and all sources (including arbitrary redshifts for remaining sources) 
[dotted line]. 
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Figure 9. Star formation history. The heavy curve is the prediction of Model S, with 
vertical lines indicating the uncertainty. The remaining lines indicate star formation 
histories predicted by several other models: solid line (Baugh et al 1998, [22]), dotted 
line (Guiderdoni et al 1999, [26]), other broken lines (Dwek et al 1998, [8], figure 3). 



